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We present a simple formalism describing evolution of a qubit in the process of its measurement 
in a circuit QED setup. When a phase-sensitive amplifier is used, the evolution depends on only 
one output quadrature, and the formalism is the same as for a broadband setup. When a phase- 
preserving amplifier is used, the qubit evolution depends on two output quadratures. In both cases 
a perfect monitoring of the qubit state and therefore a perfect quantum feedback is possible. 
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INTRODUCTION AND QUALITATIVE 
DISCUSSION 



The goal of this Lecture is to present a physical picture 
of the process of continuous quantum measurement of a 
qubit in the circuit quantum electrodynamics (cQED) 
setup [ll-Il] (Fig. 1), extending or reformulating the pre- 
vious theoretical descriptions p-Q • Understanding of the 
qubit evolution in the process of measurement is impor- 
tant for developing an intuition, which is useful in many 
cases, in particular in designing various schemes of the 
quantum feedback llCH12l|. When a quantum measure- 
ment is discussed [13|, there are usually two different 
types of questions to answer: we can either focus on ob- 
taining information on the initial state (before measure- 
ment) or focus on the quantum state after the measure- 
ment (i.e. evolution in the process of measurement). Let 
us emphasize that we consider the latter problem here 
and essentially extend the collapse postulate by describ- 
ing continuous evolution "inside" the collapse timescale. 

In the cQED setup (Fig. 1) a qubit interacts with 
a GHz-range microwave resonator, whose frequency 
slightly changes depending on whether the qubit is in the 
state |0) or |1) In turn, this frequency shift affects 

the phase (and in general amplitude) of a probing mi- 
crowave, which is transmitted through the resonator (in 
another setup the microwave is reflected from the res- 
onator, but the difference is not important). The outgo- 
ing microwave is amplified, and after that the rf signal is 
downconverted by mixing it with the original microwave 
tone, so that the low-frequency (< 100 MHz) output of 
the IQ mixer provides information on the qubit state. 
The output noise is mainly determined by the first ampli- 
fying stage, the pre-amplifier. With recent development 
of nearly quantum-limited superconducting parametric 
amplifiers [3, EBl; it is natural to use them as pre- 
amplifiers '9] instead of cryogenic high-electron-mobility 
transistors (HEMTs) [H-Ql , which usually have noise tem- 
perature above 3 K. 

Continuous quantum measurement in the cQED setup 
in some sense falls in between a qubit measurement by a 
quantum point contact (QPC) or a single-electron tran- 
sistor (SET), the theory of which has been developed over 
a decade ago [16, ITJ, and continuous quantum measure- 
ment in optics, for which the theory of quantum trajec- 
tories has been developed even earlier [1^. Neverthe- 
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FIG. 1: Schematic of the cQED setup. Microwave field of 
frequency tOm is transmitted through (or reflected from) the 
resonator of frequency uJr, which slightly changes, oj^ ± x, de- 
pending on the qubit state. After amplification the microwave 
is sent to the IQ mixer, which produces two quadrature sig- 
nals: I{t) and Q{t). For a phase-preserving amplifier we de- 
fine I{t) as the quadrature carrying information on the qubit 
state, while for a phase- sensitive amplifier we define I(t) as 
corresponding to the amplified quadrature. 



less, the cQED setup differs from both these cases, and 
this is probably the reason why there is still a confu- 
sion about the proper physical description of the mea- 
surement process. The measurement by the QPC or 
SET is of the broadband type, meaning that the mon- 
itored frequency band starts from zero. In contrast, the 
cQED setup is of the narrowband type: we deal only 
with a relatively narrow band around the probing mi- 
crowave frequency uJm- This necessarily involves two or- 
thogonal quadratures [l9j : we work with rf signals of 
the type A{t) cos{uJmt) + B[t) sin{u)mt), and there are es- 
sentially two signals A{t) and B{t) instead of only one 
in the broadband case. In this sense the cQED setup 
is similar to the optical (especially cavity QED) setup 
fisj : however, there is an important difference: in the 
cQED case the outgoing microwave is amplified (Fig. 1) 
before being mixed with the original microwave, while 
there is no amplification stage in the standard optical 
setup. The operation will obviously depend on whether 
a phase-sensitive or a phase-preserving amplifier is used, 
since a phase-preserving amplifier necessarily adds the 
half quantum of noise into any quadrature 0, [20l - [2^ . 
Notice that the quantum trajectory theory for the cQED 
setup was developed in however, the amplifier stage 
was essentially missing in the analyzed model. 

In this Lecture we consider the simplest cQED case, as- 
suming dispersive regime [6| , exactly resonant microwave 
frequency, absence of the Rabi drive, and sufficiently wide 
resonator and amplifier bandwidths for the Markov ap- 
proximation. Some generalizations are rather straight- 
forward; however, our goal is a simple picture in a simple 
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case. 

A description of continuous qubit measurement is es- 
sentially a description of the quantum back-action. Fol- 
lowing the same quantum Bayesian framework as for 
the measurement by QPC/SET ^ (see [H for re- 
view), we will discuss two kinds of measurement back- 
action onto the qubit, which we name here as "spooky" 
and "realistic". The "spooky" (or "quantum", "infor- 
mational", "non-unitary") back-action does not have a 
physical mechanism and therefore cannot be described 
by the Schrodinger equation (in contrast to what people 
often think, trying to find a mechanism for the quantum 
collapse); however, it is a common-sense consequence of 
acquiring information on the qubit state in the process 
of measurement. This is essentially the same back-action 
which is discussed in the EPR paradox jlJ] and Bell in- 
equality violation pBj : the only difference is that in our 
case the information is incomplete and therefore we have 
to use the quantum Bayes rule [H, [H, [23| instead of the 
projective collapse rule. In contrast, the "realistic" (or 
"classical" , "unitary" ) back-action has a physical mecha- 
nism: in the cQED case it is a fluctuation of the number 
of photons in the resonator, which affects the phase of 
the qubit state. The "realistic" back-action is usually 
discussed in the standard theories of the cQED measure- 
ment Actually, there is a certain spookiness even 
in the "realistic" back-action (it may be affected by a 
delayed choice, as discussed in Conclusion); however, we 
do not want to emphasize it to keep the picture simple. 
When we measure the z-coordinate of the qubit state on 
the Bloch sphere (the basis states |1) and |0) correspond 
to the North and South poles), then the "spooky" back- 
action changes the z-coordinate and leads to the state 
evolution along the meridian lines, while the "realistic" 
back-action leads to the evolution around the z-axis, i.e. 
along the parallels. 

It is important to notice that when the probing mi- 
crowave leaves the resonator after interaction with the 
qubit, one quadrature of the microwave carries informa- 
tion about the qubit state, while the orthogonal quadra- 
ture carries information on the fluctuating number of 
photons in the resonator [6 9]. Therefore, if a phase- 
preserving amplifier is used, then the "spooky" and "real- 
istic" back-actions are fully separated and correspond to 
two orthogonal quadratures I{t) and Q{t) measured after 
the mixer (it is trivial to choose the proper linear combi- 
nations of the I/Q mixer outputs). The signals I{t) and 
Q{t) are necessarily noisy, and the measurement back- 
actions are stochastic; however, there is a correlation (full 
correlation in the ideal case) between the output noise 
and the back-action noise in both channels. As a result 
(derived later), for a quantum-limited phase-preserving 
amplifier and absence of extra decoherence, the measured 
quadratures I{t) and Q{t) give us full information about 
the back-action, so that a random evolution of the qubit 
wavefunction can be monitored precisely (a useful anal- 
ogy is with a Brownian particle under a microscope: we 
cannot predict its motion, but we can monitor it). This 



is what is needed, in particular, for arranging a perfect 
quantum feedback control of the qubit state. It is inter- 
esting to notice that for an ensemble-averaged evolution 
(in which the random but monitorable qubit evolution is 
replaced by dephasing) , exactly one half of the ensemble 
dephasing F comes from the "spooky" back-action, and 
the other half comes from the "realistic" back-action. 

In the case of a phase-sensitive amplifier it is sufficient 
to measure after the mixer only the quadrature which 
was amplified; let us still denote it as I{t), though now 
its phase is determined by the amplifier instead of the 
microwave-qubit interaction. In this case the "spooky" 
and "realistic" back-actions are in general mixed (not 
separated), because there is only one output signal I{t). 
This situation exactly corresponds to the broadband 
measurement by the QPC/SET with a correlation be- 
tween the output and "realistic" back-action noises [1^ . 
The situation simplifies when the amplified quadrature is 
the one which carries information about the qubit state 
(z-coordinate). Then in the quantum-limited case the 
"realistic" back-action is fully absent: we cannot measure 
the photon number fluctuation and correspondingly it 
does not fluctuate (in the imperfect case the effect of the 
remaining "realistic" back-action can be described by an 
extra dephasing). So we are left with only the "spooky" 
back-action, and the quantum measurement description 
coincides with the simpler theory of measurement by a 
symmetric QPC [1^, which does not produce the "re- 
alistic" back-action. In contrast, in the case when the 
photon-number quadrature is amplified, we do not obtain 
any information on the qubit z-coordinate, and therefore 
there is no "spooky" back-action, but only the "realistic" 
one. In a general case, when the amplified quadrature 
makes an arbitrary angle (p with the qubit-information 
quadrature, both types of the back-action are present, 
and their strength depends on ip. It is important to 
mention that the ensemble dephasing rate F does not 
depend on ip, as required by causality. In particular, in 
the quantum-limited case the contribution F cos^ ip comes 
from the "spooky" back-action, while F sin^ tp comes from 
the "realistic" back-action. 

Let us emphasize that both the phase-sensitive and 
phase-preserving amplifiers permit exact monitoring of 
the qubit state and therefore a perfect quantum feed- 
back. The necessary condition in both cases is that the 
detection system is quantum-limited. 

In the following sections a formal description of the 
above discussed results is presented. We start with re- 
viewing the Bayesian approach for the broadband qubit 
measurement, then briefiy discuss the difference between 
phase-preserving and phase-sensitive amplifiers, and then 
present the formalism of the narrowband continuous mea- 
surement of a qubit in the cQED setup. In Conclusion we 
briefly discuss generalizations of the formalism, quantum 
feedback, and the causality principle. We note that our 
approach can be converted into the formal language of 
the positive-operator-valued-measure(POVM)-type gen- 
eralized quantum measurement [28| (then separation of 
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FIG. 2: Schematic of a broadband measurement setup: 
double-quantum-dot qubit is measured by a QPC (tunnel 
junction). The output signal I{t) is the QPC current. 



the "spooky" and "realistic" back-actions corresponds to 
the decomposition of the measurement operator into di- 
agonal and unitary parts - see later) , and our results for 
the case of a phase-sensitive amplifier are very similar to 
the results of Ref. @ . 



II. BROADBAND MEASUREMENT 

In this section we review the Bayesian formalism 
[l^H^ for the broadband measurement of a qubit, con- 
sidering only the simple case without additional evolu- 
tion, and thus emphasizing the main physical idea of the 
formalism. We start with the broadband formalism be- 
cause it is simpler than for the narrowband (cQED) mea- 
surement and it can be used as a natural step in under- 
standing the cQED setup. 

For definiteness let us assume that the qubit is a dou- 
ble quantum dot populated with one electron (Fig. 2), 
and the states |0) and |1) correspond to the electron lo- 
calized in one or the other dot. The qubit is measured 
by a small-transparency tunnel junction (model of QPC), 
whose barrier height depends on the electron location, so 
that the two qubit states correspond to different average 
currents Iq and Ii through the QPC. The voltage across 
the QPC is sufficiently large to make the detector out- 
put classical (Markov approximation), and |A/| ^ |/c|, 
where AI = Ii — Iq is the response and Ic = {lo + -^i)/2 
is the mean value; this weak response assumption allows 
us to consider the QPC current I{t) as a quasicontinu- 
ous noisy signal (see [23,] for the detailed discussion of 
required assumptions; the formalism needs only a minor 
change if A/ ^ Ic)- Then the output signal of the detec- 
tor is 

I{t)=h + {M/2)z{t)+^{t), S^{oj)^S, (1) 

where z = pn — poo is the z-component of the Bloch 
sphere representation of the qubit density matrix p{t), 
and ^{t) is the white shot noise with spectral density 
S = 2elc (we use the single-sided definition for the spec- 
tral density, in which the signal variance ("power") cor- 
responds to S'(w) dw/27r; the definition of S is twice 
smaller in Ref. [?! and Att times smaller in Refs. [8l[27j). 
We emphasize that the detector signal /(t) is classical, 
and the qubit state p{t) is practically unentangled from 
the detector, but obviously depends on I{t). 



The detector Hamiltonian and the qubit-detector in- 
teraction Hamiltonian are given in Refs. [l^, [2^, they 
are not really important for our discussion here. For 
simplicity let us assume that the qubit Hamiltonian is 
zero, Hqb = 0, so that the qubit evolution is due to the 
measurement only. In this case the qubit evolution dur- 
ing time t happens to be determined only by the time- 
averaged value of the measured detector output 



Irnit) = -[ I{t')dt', 
t Jo 



(2) 



which would contain full information for a classical mea- 
surement. 

Because of the correspondence principle, the evolution 
of the diagonal elements of the qubit density matrix p 
(Trp — 1) should correspond to the classical evolution 
of probabilities, which are given by the classical Bayes 
rule. The Bayes rule says that an updated (a posteriori) 
probability of a system state is proportional to the initial 
(a priori) probability and the probability (likelihood) of 
the obtained measurement result assuming this particu- 
lar state. In our case Im{i) is the measurement result, 
and its probability for the qubit in the basis state | j) has 
the Gaussian distribution 

P|,)(7„) = -jL= exph(7„ - Ijf/2D], Z? - |, (3) 

where D is the variance, which decreases with the mea- 
surement time t. Therefore the correspondence principle 
demands the Bayesian evolution 



Piiit) _ Pii(O) exp[-(/™(t)-7i)V2i?] 



Poo(t) Poo(O) exp[-(/,„(0 - hfl2D] ' 



(4) 



which in our terminology is due to the "spooky" back- 
action; it cannot be described by the Schrodinger equa- 
tion, but follows from common sense. 

If the phase of qubit state is not affected by 
the measurement process (no "realistic" back- 
action), then an arbitrary initial wavefunction 
IV'(O)) = VM O)|0) + e^VmMll) becomes 

= Vpooit) |0) + e^^VPiiWIl) with the same 
phase 0; therefore for an arbitrary mixed state we get 



Poi{t) = Poi(O) 



\^ Poojt) piijt) 
Vpoo(0)pii(0) 



(5) 



Equations (|4]) and (O describe the "spooky" back-action. 

Now assume that due to the qubit-detector interaction 
(e.g. Coulomb interaction), each electron passing through 
the detector rotates the qubit phase 4>hy a small amount 
A(f). From the measured result Im{t) we know exactly 
how many electrons passed through [n^ = Im{t)t/q with 
q being the electron charge] , and can easily introduce the 
corresponding phase factor into Eq. ([S]) : 



Poi {t} = Poi (0) = exp[iKI,n (t) t], (6) 

\/Poo(0)pii(0) 
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where K = A(j>/q. The non-stochastic factor exp{iKIct) 
can be obviously ascribed to the qubit Hamiltonian; 
however, this is not important here. The factor 
ex-p[iKIm{t) t] in Eq. ^ is the effect of the "reahstic" 
back-action. It may or may not be present in a particular 
physical situation; for example, K ~ ioT measurement 
by a symmetric QPC, while 7^ in an asymmetric 
QPC or SET case. 

Finally, if there is an extra pure dephasing of a qubit 
with rate 7, then Eq. © becomes 



Poi{t) = PaiiO) ^ == exp[iKI^{t)t]e ^ . (7) 

VPoo(O) Pii(O) 

Equations (U) and ([7]) is the main starting point of 
the Bayesian formalism [23| . It is then easy to in- 
clude non-zero qubit Hamiltonian Hqb by differentiat- 
ing Eqs. (H]) and ([7]) over time (paying attention to 
whether the Stratonovich or Ito definition of the deriva- 
tive is used) and adding terms due to i/qj,. Energy relax- 
ation and other mechanisms of the qubit state evolution 
can be included in the same way. Actually, there are 
many ways to derive the Bayesian equations @ and ([7]) 
[3 US [Hj H^l , but we focus here only on their mean- 
ing, not on their derivation. 

Notice that averaging of Eqs. (U) and (O over the mea- 
surement result Im (i.e. ensemble averaging) with the 
probability distribution 



P{Irn) = poo(0)P|o)(/m) +Pll(0)P|l)(/. 



(8) 



gives the same evolution as for a pure dephasing: the 
diagonal matrix elements of p do not evolve, while the 
off-diagonal element poi decays as poi(0)e~'"* with the 
ensemble dephasing rate [1^ 



(All 
45 



K^S 



+ 7, 



(9) 



which has clear contributions from the "spooky" back- 
action, "realistic" back-action, and additional dephasing. 

In the case 7 = an initially pure qubit state re- 
mains pure; in other words we can monitor evolution of 
a qubit wavefunction. This property can be used as the 
definition of a quantum-limited detector (T^, The 
quantum efficiency rj can then be naturally defined as 
77 = 1 — 7/r. If by some reason the "realistic" back- 
action is considered as dephasing (i.e. only in the aver- 
aged way), then the quantum efficiency can be defined 
as 77 = 1 — 7/r — K^S/AT (here the definitions of fj and 
rj are exchanged compared with the definitions in 23]). 
In other words, fj = (A/)^/4S'r is the relative contribu- 
tion of only the "spooky" back-action in the ensemble 
dephasing F. In particular, this definition is relevant to 
the peak-to-pedestal ratio of the Rabi spectral peak [3l| , 
which is equal to Ai]. As an example, if 7 = and con- 
tributions in Eq. ([9|) from the "spooky" and "realistic" 
back-actions are equal to each other (as in the cQED 
setup with a phase-preserving amplifier) , then 77 = 1 but 
7? =1/2. 



A non-ideal detector (77 < 1) can be modeled in two 
equivalent ways [s^ : we either add an extra dephasing 
7 to the qubit or we add an extra noise to the output 
of the ideal detector. Only the total dephasing F, re- 
sponse A/, total output noise S, and correlation factor 
K = 5{4)) /5{Imt) are the physical (i.e. experimentally 
measurable) parameters, while distribution of the non- 
ideality between the extra dephasing and additional out- 
put noise is a matter of convenience [here (j) = a'rg(poi); 
and notation (0) reminds about averaging over additional 
classical noise at the output]. We emphasize that the 
Bayesian formalism deals only with the experimentally 
measurable parameters AJ, S*, K, F, and the output sig- 
nal I{t). 

In the ideal case (7; = 1) the evolution equations ^ 
and ([5]) can be translated into the language of POVM- 
type generalized measurement. In this approach the ef- 
fect of measurement is described as [1^ 

l^(^)>=PWW ^^*)^Tr[AftM,p(0)]' ^''^ 

where Mr is the so-called measurement (Kraus) oper- 
ator, corresponding to the result R. The probability of 
the result i? is Pfl = ||A/i?|?/;(0))|p using wavefunctions or 
Pr = Tr[M|j,Mflp(0)] using density matrices; therefore 
the POVM elements Mj^M^ should satisfy the complete- 
ness condition J2r^r-^R — 1- The relation between 
this approach and the quantum Bayesian approach can 
be understood via the operator decomposition 



Mr = Ur^mIMr, 



(11) 



where Ur \s unitary and the square root of the posi- 
tive operator M^^Mr is defined in the natural way in 

the diagonalizing basis. It is easy to see that \J M^^Mr 

is essentially the quantum Bayes rule (in the diagonal- 
izing basis); in our terminology it corresponds to the 
"spooky" back-action, while Ur corresponds to the "re- 
alistic" back-action. For the discussed setup the result R 
is I rait), the "spooky" back-action [m\^MrY/'^ should be 
determined by the probabilities P\j) {Im) given by Eq. ([3|), 
and the "realistic" back-action Ur is given by the phase 
factor in Eq. (jS]). Therefore the corresponding measure- 
ment operator is 



M(/„) = exp{-iKI^taj2) 



^1o)(/m)|0)(0| + JP|i)(/™)|l)(l 



[12) 



where a, is the Pauli matrix. 



III. PHASE-PRESERVING VS. 
PHASE-SENSITIVE AMPLIFIERS 

Before discussing microwave amplifiers, let us consider 
a measurement of an oscillator, for example, a mechan- 
ical resonator with frequency oj,- and mass 772. This is a 
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very well studied problem 0, [T^l : so we will only discuss 
a way to understand the results. A classical resonator 
position X oscillates as Xc{t) = Acos{uJrt) + B8iTi{ujrt) 
(in this section x stands for the usual spatial coordinate, 
not for the Bloch sphere coordinate) . The corresponding 
quantum state is called the "coherent state" in the optical 
language; it is represented by the wavefunction ipix, t) = 
ipgr[x — Xc{t)]eicp{ipcx/h), where ipgrix) is the ground 
state and Pc — mxc{t) is the classical momentum. So the 
coherent state is essentially the ground state with oscil- 
lating center position. Notice that continuous quantum 
measurement of a resonator position can be described 
in the same Bayesian way [s^ as in the previous section; 
for example the "spooky" back-action gives the evolution 
ip{x,t) = ■!/'(a;,0)exp[-(T™(i)-J(a;))V4£']/Norm, where 
I{x) is the average detector signal for the resonator po- 
sition x, and Norm is normalization [see Eqs. ((4]) and 
([5])]. The time step t in this case should be chosen much 
shorter than ui~^ so that the unitary evolution and evo- 
lution due to measurement may be simply added. 

Let us consider the following game: Charlie prepares 
an oscillator in a coherent state with quadratures A and 
B, gives it to David, and David's goal is to find A as 
accurately as possible. An optimal strategy is rather ob- 
vious: David should make a projective measurement of 
x at time t = 2TTn/uJr with any integer n (to avoid con- 
tribution from the B-term), and the measurement result 
is the best estimate of A [if the measurement is done at 
t = (27m -|- T:)/LUr, then the result should be multiplied 
by —1]. Even though the strategy is optimal, the inac- 
curacy of David's result is obviously the width (standard 
deviation) agr = \/ h/2mijjr of the ground state shape 
\ipgr{x)\'^] in energy units this inaccuracy corresponds to 
one half of the energy quantum. 

Now assume that David cannot make projective mea- 
surements, but only "finite-strength" (i.e. imprecise) 
measurements. The best accuracy agr can still be 
achieved if the measurement is done in the simple 
but very clever "quantum non-demolition" (QND) way: 
many finite-strength measurements are made at times 
t = 2Tm/ijJr] this is called "stroboscopic" measurement 
[isj . Since the oscillator returns to the same state after 
the period 27r/wr, the unitary evolution is not impor- 
tant, and many finite-strength measurements (described 
by the Bayesian equation above) are "stacked" to pro- 
duce a strong, essentially projective measurement. More 
generally, the necessary condition to have the best accu- 
racy Ggr for A is that the measurement is not sensitive 
to the quadrature B. 

Now assume that David is only allowed to make a con- 
tinuous measurement with unmodulated weak strength 
(so that the inaccuracy achieved after io~^ is much larger 
than agr). Then the "spooky" back-action gets mixed 
with the unitary evolution, essentially adding noise into 
the monitored evolution, so that after a while the res- 
onator state becomes mostly determined by the back- 
action and almost not dependent on the initial state. As 
the result, the best accuracy for measurement of A be- 



comes yplGgri which in energy units corresponds to two 
half-quanta [13|J, that is twice worse than for the projec- 
tive or stroboscopic measurement. However, the contin- 
uous monitoring gives us an information about B in the 
same way as for A, so the accuracy of i?-measurement is 
also yplugr- Therefore, in some sense continuous phase- 
insensitive measurement brings the same total informa- 
tion as the phase-sensitive (e.g. stroboscopic) measure- 
ment; however, in our game only half of this information 
is useful for the David's goal. 

After discussing measurement of the resonator state it 
is easy to understand the quantum limits for the high- 
gain microwave amplifiers. Now suppose Charlie pre- 
pares a coherent state of a microwave resonator with 
quadratures A and B, gives it to David to find A^ 
and David uses an amplifier for amplification of the mi- 
crowave field, which slowly leaks from the resonator until 
it is empty. There is only classical signal processing af- 
ter the amplifier, so amplification is essentially the quan- 
tum measurement. The results are the same as above 
0, [13, [m, [l^l: a phase-sensitive amplifier, which ampli- 
fies only A-quadrature and "de-amplifies" (attenuates) 
-B-quadrature can measure A with accuracy Ogr, while a 
phase-preserving amplifier can measure A only with ac- 
curacy cTgr (and also measures B with the same accu- 
racy \f2 (jgr ) . Technically, the accuracy is limited by the 
noise at the amplifier output, so this noise should forbid 
measuring of A with accuracy better than agr by a phase- 
sensitive amplifier, and better than \f2agr by a phase- 
preserving amplifier. Therefore in the quantum-limited 
case the output noise power of a phase-preserving am- 
plifier (per quadrature) is twice larger than for a phase- 
sensitive amplifier with the same gain; this is often called 
an "additional noise" , corresponding to one half of the en- 
ergy quantum (one more half-quantum is present in both 
cases) 0, El, El, [H- It may be somewhat confusing 
why this result does not depend on the rate with which 
the microwave leaks from the resonator. So let us check 
the scaling: for fc times slower leakage (A; times larger 
Q-factor) the microwave amplitude is ^fk times smaller, 
but accumulation time is fc times longer, therefore the 
measured signal for the quadrature A is \fk times larger, 
which is the same factor as for the noise accumulation. 
Therefore, the signal-to-noise ratio which determines the 
A-accuracy does not depend on the leakage rate (res- 
onator bandwidth). 

IV. NARROWBAND (CQED) MEASUREMENT 

Using the above discussion of microwave amplifiers, it 
is easy to extend the Bayesian approach for a broadband 
quantum measurement to the narrowband cQED setup. 

We consider the standard cQED setup [1-9], in which 
a qubit interacts with a microwave resonator, and assume 
the dispersive regime with the Hamiltonian 

H — {hjjq}j/2) + hi/JrO^ a + hxa^ aoz^ (13) 
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where ujqi, = ujyb,bare + X is the Lamb-shifted qubit 
frequency with no photons in the resonator, x — 
/ (ujqbfiare — w,) is the effective couphng with g be- 
ing the Jaynes-Cummings couphng, uir is the bare res- 
onator frequency, Pauh operator Uz acts on the qubit 
state in the energy basis {|0), and the resonator cre- 
ation/annihilation operators are at and a. Notice that 
the resonator frequency increases by 2% when the qubit 
state changes from |0) to |1); conversely, the qubit fre- 
quency increases by 2% per each additional photon in the 
resonator. To measure the qubit state, a microwave field 
with frequency Wm is either transmitted through or re- 
flected from the resonator, then amplified and sent to the 
IQ mixer, which measures both quadratures relative to 
the original microwave tone (Fig. 1). The qubit state af- 
fects the resonator frequency and therefore affects the 
phase (and in general amplitude) of the outgoing mi- 
crowave. 

An elementary Fabry-Perot analysis gives the classical 
(complex) microwave field Fj. inside the resonator: 

= Z TT7 r-j-, 

1 — ZiyuJm — (^r)l 1^ 

where is the applied incident field, tin is the trans- 
mission amplitude of the barrier from the incident side, 
K is the resonator bandwidth due to the microwave leak- 
age from the both sides (the Q-factor is w^/k), and the 
round-trip time is r^t — 2'K/uJr for a half- wavelength 
resonator and r^t — tt/w^ for a quarter- wavelength res- 
onator. A similar formula with the same denominator 
describes a lumped resonator. In presence of the qubit, 
the resonator frequency lo^ in this formula is substituted 
by bjj. ± X, depending on the qubit state. Notice that 
for the quantum measurement analysis there is no dif- 
ference between the cases of transmission and reflection 
for the same Fr and k, because the field leaking from 
the resonator is determined only by Fr and k. (The re- 
flection case has a technical advantage of dealing with a 
twice smaller outgoing microwave field for the same mea- 
sured signal.) However, an important parameter is the 
collected fraction r]coi = Kcoi/k of the leaking microwave 
power; we will often assume the ideal case r]coi = 1 (for 
the transmission setup this requires strongly asymmetric 
coupling, |t„J < \tout\)- 

For simplicity we assume the resonant case, ajm — 
then the ensemble qubit dephasing due to measurement 
is in 

r = 8x'n/K, (15) 

where n is the average number of photons in the res- 
onator. It is easy to include Rabi oscillations into the 
model; however, we do not do it for simplicity and also 
for more transparent analogy with Sec. |TT1 in which we 
considered a qubit with zero Hamiltonian, evolving only 
due to measurement; this case exactly corresponds to the 
cQED Hamiltonian ([T^ in the rotating frame. 

We will need several assumptions to describe the qubit 
state evolution in the process of measurement. First, for 



the validity of the dispersive approximation (jl3p we need 
sufficiently large qubit- resonator detuning, jw^fc — w^l ^ 
and not too many photons in the resonator, n <^ 
{u!qb—ujr)'^/g'^ (wc do not consider the recently discovered 
nonlinear regime [13]). Second, to use the Markov ap- 
proximation for the evolution we need the so-called "bad 
cavity" assumption: F <C k ^ (if the qubit evolves 
due to Rabi oscillations with frequency VIr, we also need 
^ k). This assumption means that the photons leave 
the resonator much faster than evolution of the qubit 
state, and therefore there is practically no entanglement 
between the qubit and unmeasured microwave field. This 
assumption also implies that the two resonator states for 
the qubit states |0) and |1) are almost indistinguishable, 
h{x/ k)'^ <C 1. Third, we use the "weak response" as- 
sumption, which requires a small phase difference be- 
tween the two resonator states, |xl/'* ^ 1- This means 
that each outgoing photon carries only a little informa- 
tion about the qubit state. Notice that for n > 1 the 
previous assumption k ^ T automatically implies the 
weak response, and even for n <C 1 the weak response 
assumption is not always needed. Fourth, we will neglect 
the qubit energy relaxation due to measurement [g, 0] , 
which can be added later. 

A coherent state in the resonator with average n 
photons and zero average phase corresponds to the 
oscillation of the field expectation value {Fr{t)) — 
cos{uJmt), where agr is the ground state width 
(rms uncertainty) and we assume w„j = ujr- (Notice 
that the amplitude agr corresponds to 1/4 photon.) 
Interaction with the qubit slightly changes the phase, 
cos{u}mt T Sx/'^), depending on the qubit state, so that 

{Fr{t)) = Acos{ujmt) + B sm{ujmt), A — 2\ffiagr, 
B = ±(4x/K)%/^crgr = (4x/K)Vna-gr z, (16) 

where z is the qubit Bloch coordinate. Thus the small 
_B-quadrature carries information about the qubit state, 
while larger A-quadrature may give us information on the 
fiuctuations of the photon number in the resonator. In 
the optical representation (Fig. 3) with axes A/ugr and 
Bjogr, the two resonator states for the qubit states |0) 
and |1) are shown as two "error circles" 35] with rms un- 
certainty 1 along any direction and distance 2-\/H between 
the origin and circle centers (if axes Ajlogr and B jlogr 
are used, then the distance is -y/ffi while the uncertainty 
is 1/2). 

A. Phase-sensitive amplifier 

Let us start with the case when a phase-sensitive am- 
plifier is used in the cQED setup. Also, we first assume 
the most ideal case: the amplifier is quantum-limited, 
it amplifies the optimal _B-quadrature, there is no mi- 
crowave collection loss {kcoi/h = 1), and there is no ex- 
tra noise or dephasing. Then, as discussed in the pre- 
vious section, measuring once the microwave contents of 
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FIG. 3: Phase space representation: for each qubit state the 
coherent state with (Fr) — Acos{ijJrt) + Bsm{iUrt) in the res- 
onator is shown [sH] as an "error circle" with radius 1, shifted 
by from the origin. Axes are normalized by the stan- 

dard deviation agr of the ground state. The _B-quadrature 
carries information about the qubit state, the ^-quadrature 
corresponds to the number of photons in the resonator. 




phase-sensitive (a) 




phase-preserving (b) 



FIG. 4: Relations between relevant quadratures. The quadra- 
tures A and B are for the microwave field in the resonator, 
(a) For a phase-sensitive amplifier the amplified quadrature 
makes an angle ip with the informational B-quadrature. The 
corresponding quadrature at the mixer output is defined as 
I{t) [the output Q{t) is then useless]. (b) For a phase- 
preserving amplifier we define output quadratures I{t) and 
Q{t) as corresponding to the resonator quadratures B and A. 



the resonator (by fully emptying it) we can measure the 
-B-quadrature with imprecision agr- Therefore, in the 
continuous measurement for time t the B-quadrature is 
measured with imprecision o-gr/v^, which converts into 
the imprecision K/t/{AxV^) of the qubit z-coordinate. 
Following the language of Sec. HIl let us discuss the signal 
and noise at the output of the setup. There are two out- 
puts of the IQ mixer; however, only the amplified quadra- 
ture carries an information, so let us denote the corre- 
sponding output of the mixer (or their linear combina- 
tion) as I{t). Then the response AI = Ii — Iq corresponds 
to Az — 2 and AB = 8{x/ K)V^<^gr- For measurement 
during time t the above variance {k / 1) / {Ax-\/n)^ of the z- 
coordinate converts into the variance (K;/t)(A//8xV^)^ 
of the measured output Im = (1/*) JoI{t')dt'. Equat- 
ing it with D = S/2t, we find the (single-sided) spectral 
density of the I{t) noise: 



(AW)2K/(32x'n), 



(17) 



where we replaced S with Smin and AI with AI„iax to re- 
mind that we consider the quantum-limited case, and the 
response is maximized by amplifying the optimal quadra- 
ture. Notice that since Almax oc the noise Smin 
does not depend on the qubit or resonator properties; 
it is essentially the amplified vacuum noise and depends 
only on the amplifier gain. 



Obtaining information on the qubit z-coordinate via 
the signal I{t) with response Al^ax and noise Smin, we 
necessarily cause the "spooky" back-action described by 
Eqs. @ and ([S]). As discussed in Sec. [TTl this is the 
consequence of the corresponding principle or just the 
common sense. Now averaging the poi evolution in Eq. 
Q over the measurement result /,„ with its probability 
distribution ([5]) and we see that the "spooky" back- 
action dephases an ensemble of qubits with the rate [see 
Eq. ^] {AlmaxY / ^Smin = Sx^^/^. This rate coincides 
with the total ensemble dephasing (|15p . and therefore 
the qubit state cannot additionally fluctuate due to any 
other reason. Thus we derived an important result: in the 
ideal case with phase-sensitive amplifier there is only the 
"spooky" hack- action and no "realistic" hack- action. This 
means that the numher of photons in the resonator does 
not fluctuate (otherwise there would be an additional de- 
phasing), which makes sense since we cannot measure 
the A-quadrature, carrying information on the photon 
number. Notice that it is also easy to prove this result 
when ujm ^r- Then from Eq. (|14p we obtain that the 
informational quadrature amplitude is multiplied by the 
factor [l-|-4(w,„ — '^r)^/^^]"^^^ compared with Eq. 
The response Almax is multiplied by same factor, while 
the noise Smin does not change. Therefore, the "spooky" 
back-action contribution into the ensemble dephasing is 
multiplied by the factor [1 -I- Aiujm — '^r)^/^^]"^, which 
again coincides with the result [1, 0| for the total ensem- 
ble dephasing F. This proves the absence of the "real- 
istic" back-action for the non-resonant case uim 7^ as 
well. 

Now let us consider the case when an ideal phase- 
sensitive amplifier amplifies the A-quadrature (we again 
assume ujm — for simplicity). Then we do not get 
any information on the qubit z-coordinate, and therefore 
there is no "spooky" back-action, but there is the "real- 
istic" back-action due to fluctuating number of photons. 
The description of evolution in this case is essentially the 
standard description [1, 0]- Let us still denote with I{t) 
the output signal from the mixer, corresponding to the 
amplified quadrature. For measurement during time t we 
measure A-quadrature with imprecision agr/^/Kt. This 
is consistent with fiuctuation of the number TV of emit- 
ted photons: var(7V) = N, N = nnt. The correlation 
function of photon number in the resonator depends on 
time as exp(— Ki/2) @, 0|) which means that each extra 
photon inferred from I(t) fluctuation spends (on aver- 
age) time 2/k in the resonator and therefore changes the 
qubit phase (j) by Ax/k (the correlation time 2/k is es- 
sentially the lifetime of the field, not power @, Q)- Then 
0-variance is var(0) = [Ax/ K^nKt, and the correspond- 
ing ensemble dephasing is va.r{(j>)/2t — %x^njK. As ex- 
pected, this reproduces the standard result ([T5|) for the 
ensemble dephasing, while for individual qubit evolution 
we have the above discussed correlation: each additional 
photon inferred from /(i) fluctuation changes by 4x/k. 
For the same amplifier gain and noise as for measuring 
B-quadrature, we get = (AxJ n)^{blm) I Almax, 
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and therefore the correlation is K = S{(f>) /d{Imt) = 
32{x^/K)n/ Almax = Almax/Smin- It is easy to check 
that K^Smin/'i [see Eq. ^] coincides with ensemble de- 
phasing r from Eq. ([TS]), as expected for the presence of 
only "realistic" back-action. 

Finally, assume that the phase-sensitive amplifier am- 
plifies the quadrature, which makes angle ip with the 
optimal _B-quadrature and angle 7r/2 — ip with the A- 
quadrature. The measured signal I{t) still denotes the 
output of the IQ mixer, corresponding to the amplified 
quadrature; now it gives information about both B and 
A quadratures, with the factors cos f and sin ip, respec- 
tively. Combining the "spooky" and "realistic" back- 
actions, we get the same formulas as for the broadband 
detection of Sec. |lTl 



piiW pii(o) 



exp[/„(t)A//i?], D 



2t' 



eyi^[iKim{t) t] 



(18) 
(19) 



Pm(t) poo(O) 
Pm{t) _ Vpoajt) pii{t) 
Poi(O) " v/poo(0)pii(0) 

A/ = /i - /o = A/maa; COS ifi, K = K^ax siu Lp, (20) 
1 J A- T 

Ut)^- I{t')dt'-^^, (22) 

Ht) = + ^ m + m, s^ic) ^ s. (23) 

Here we introduced Im by subtracting the constant (Iq + 
/i)/2 from /,„ and did a simple algebra to convert Eq. 
(jl]) into Eq. (fT8|) : the qubit rotating frame corresponds to 
n photons, Kmax is the above discussed correlation for 
v4-quadrature amplification, \Imax is the response for 
-B-quadrature amplification, and S = Smin- Notice that 
the total ensemble dephasing ^ does not depend on ip: 

{Aljnax cos ip)^/AS„iin+ [Kmax sin ip)'^Srmn/ 4. = E. (24) 

So far we discussed only the ideal case. There are sev- 
eral mechanisms for non-ideality. First, the qubit may 
have additional environmental dephasing jenv This will 
lead to the extra factor e~'''=""' in Eq. and increase 
the ensemble dephasing F by "fenv Following the defini- 
tions in Sec. |lTl the corresponding quantum efficiency is 
Venv — (1 + JenvK/8x'^n)~^ . Second, not all microwave 
power leaking from the resonator may be collected and 
amplified. This can be characterized by the collection 
efficiency rjcoi — Kcoi/k and multiplies the response A/ 
and correlation K by the factor ^/rjcoi , while not affecting 
the output noise S. Third, if the phase-preserving am- 
plifier is not quantum-limited, it introduces additional 
noise Sadd compared with the quantum limit Smin [given 
by Eq. (|17|) when tJcoI = 1]- The corresponding ampli- 
fier efficiency is riamp = Smin/ {Smin + Sadd)- This does 
not affect A/ but multiplies K by rjamp (because for un- 
correlated Gaussian-distributed random numbers xi and 
X2^ the averaging of xi for a fixed sum xi -\- X2 gives the 
correlation (a;i)/(xi -|-a;2) = var(a;i)/[var(a;i) -|-var(a;2)]). 



If all three mechanisms of the non-ideality are present, 
then the evolution can still be described [s^] by Eqs. P^ - 
(|23l) , but 5* is now the total (experimental) output noise, 
Almax is the experimental response for = 0, so that 
S = {AImax)'^n/ (32x^11- r]coir]amp), the correlation K = 
5{4>)/6{imt) is still given by Eqs. (HOll-dlll), and the only 
change is the extra factor in Eq. (fT9|) : 



- VPoo(<)Piift)_ ,,^u j<T (t)t]e' 

Poi(O) Vpoo(0)pii(0) 

7 = F - {AImax)V4.S, 



-It 



(25) 
(26) 



where the ensemble dephasing is now F — 8x^^/'« + 7enD- 
We emphasize that the qubit evolution depends only on 
the experimentally measurable parameters Almax, S, F, 
ip, and the output signal I{t). 

The quantum efficiencies can be expressed as 



V = 



AST 



Vamp Vcol Venvi V ^ V COS (p, 



(27) 



where as in Sec.|TIl 77 is the relative contribution to F from 
both the "spooky" and "realistic" back-actions, while 
77 is the relative contribution from only the "spooky" 
back-action. The definition of fj corresponds to re- 
placing the "realistic" back-action factor ex-p{iKImt) 
in Eq. (j25p with the corresponding ensemble dephasing 
exp{—K'^St/4:). As mentioned in Sec. |lTl the peak-to- 
pedestal ratio of the spectral peak of contiuous Rabi os- 
cillations is 4fj = 4?7COS^ (p. 



B. Phase-preserving amplifier 

Now assume that a phase-preserving amplifier is used 
(this includes parametric amplifier, HEMT, etc.). Now 
both the A-quadrature and B-quadrature of Eq. 
are amplified independently with the same gain. Cor- 
respondingly, both quadratures at the IQ mixer output 
carry physical information instead of only one quadra- 
ture in the case of a phase-sensitive amplifier. Let us 
denote with I{t) the output of the IQ mixer, correspond- 
ing to the i?-quadrature; thus I{t) provides information 
on the qubit z-coordinate. The output signal for the or- 
thogonal quadrature is denoted Q{t); it corresponds to 
the A-quadrature in the resonator and provides informa- 
tion on the fluctuating number of photons. The main 
difference from the case of a phase-sensitive amplifier is 
that now the "spooky" and "realistic" back-actions are 
related to two different output signals: I{t) and Q{t). 

Let us start with the quantum-limited case and assume 
an amplifier with the same gain as in the phase-sensitive 
case, so that the /(t)-channel response is the same as 
the optimal phase-sensitive response, AI = Almax- The 
"spooky" back-action is always described by the quantum 
Bayes formulas (H))-®, but now the noise S of the out- 
put I{t) is twice larger than the value (IT71) for the phase- 
sensitive amplifier, S = 2Smin (see discussion in Sec. lIIII) : 
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therefore the "spooky" evolution is twice slower than in 
the phase-sensitive case with ip — 0. The signal Q{t) 
has the same noise S = 2Smin', it is again twice larger 
than for the phase-sensitive case with tp = 7r/2; therefore 
the correlation factor K = S{(j)) /S[J* Q{t') dt'] for the 
"realistic" back-action is twice smaller: K — Kmax/'2 
(this reduction is similar to the effect of a non-ideal am- 
plifier discussed above). We see that K — AI/S, and 
the ensemble dephasing is at least (A/)^/4S' + K^S/4 — 
(A/)^/25 = {AIjnax)/4:Smin- This again coincides with 
r = Sx^n/n, and therefore there can be no additional 
evolution of the qubit besides these "spooky" and "real- 
istic" back-actions. 

Thus in the ideal case the qubit evolution is 

'-^,-^,eMLmi/D], D = l- (28) 
Poo(*) Poo(O) 2t 



Poi{t) _ VPoo{t) pii{t) 



Poi(O) Vpoo(0)pii(0) 
A/ = /i-/o, K = AI/S, 

- ,(t) = - / Qit')dt' ~{Q), Sq 



t 



exp[iXg„(t)t], (29) 
(30) 

Si = S*, (31) 



where {Q) is the average value of Q{t) (which depends 
on n), Im{t) is defined by Eq. (|22|) . and the channels 
I{t) and Q{t) both have the same (uncorrected) noise 
S = (A/)2«;/(16x^n). Notice that (A7)V4S' = K'^S/A = 
'^X^n/n, and therefore in the phase-preserving case the 
ensemble dephasing T contains equal contributions T/2 
from the "spooky" and "realistic" back-actions. We em- 
phasize that Eqs. (j28p ~ (|29p still allow us to monitor a 
qubit wavefunction if the initial qubit state is pure. 

A non-ideal case can be analyzed in the same way as for 
the phase-sensitive amplifier. An extra dephasing "fenv of 



the qubit is described by rje 



(l + 7e™K/8x n)" 



perfect collection efficiency is described by rjcoi = Kcoi/k, 
and the amplifier efficiency is rjamp- We define rjamp — 
Sqi/S for a phase-preserving amplifier by comparing its 
output noise S (per quadrature) with the quantum limit 
for a phase-preserving amplifier: 5*^; = 2Smim so that 
Vamp = 1 in the quantum-limited case. We also define 
f] — Smin/S by comparison with a phase-sensitive am- 
plifier having the same gain, so that rjamp = Vamp/2 and 
obviously fjamp < 1/2. Similarly to the phase-preserving 
case, incomplete microwave collection multiplies the re- 
sponse AI and correlation K by the factor ^Jrjcoi but does 
not change the noise S] the extra noise in the amplifier 
multiplies K by rjamp but does not change AI. 

The qubit evolution can still be described by Eqs. 
([3T|) with the only change in Eq. p9|: 



Poi(^) ^ \Jpoa{t) Piijt) 
Poi(O) ' Vpoo(0)pii(0) 
7 r - 2(A/)V45, 



eMiKQm{t)t]e-<\ (32) 
(33) 



jenv is the total ensemble dephasing. The qubit evolution 
is determined by the parameters AI, S, T, and output 
signals I{t) and Q{t). 

The quantum efficiencies are 

V = VampVcoiVenv = (A/)^/25r, fj = rj/2. (34) 

Here rj is the fraction of F due to the contribution from 
both the "spooky" and "realistic" back-actions. The ef- 
ficiency 77 = VampVcoiVenv is the fractiou from only the 
"spooky" contribution; it corresponds to replacing the 
term exp{iKQ„i t) in Eq. p2p with the dephasing term 
exp(— ii'^t/45'). In particular, the peak-to-pedestal ratio 
of the Rabi spectral peak for the signal I{t) is Afj — 2ri. 

Let us mention that Eqs. ([28l) -([30 |) for the ideal phase- 
preserving case can also be obtained from Eqs. (fT8|)- 
(l?T|) for the phase-sensitive case in the following way. 
Let us think about a phase-preserving amplifier as a 
phase-sensitive amplifier, in which the angle ip rapidly 
changes with time, and we have to average over ip. 
When coefficients cosi^ and simp in Eq. (|20l) are sub- 
stituted into Eqs. (fT8|l and ([T9|. we see a natural for- 
mation of the quadratures Im and Qm of the phase- 
preserving setup. Then the exponential factor in Eq. 
([T^ becomes exp{ImAImax/D), and the exponential fac- 
tor in Eq. p9)) becomes exp{iKmaxQmt). Now let us 
take i nto account that the average response is AI = 
cos'^ ip Almax — AImax/2, and the phase-sensitive am- 
plifier noise S splits equally between the I(t) and Q(t) 
quadratures (the orthogonal, de-amplified quadrature is 
noiseless). The mutual cancellation of these two factors 
of 2 leads to the same form of Eq. ([25]) as in Eq. (|18p and 
the relation K = AI / S in Eq. (150)) . 

One more way to understand the relation between the 
ideal phase-sensitive and phase-preserving cases is the 
following. Instead of using a phase-preserving amplifier, 
let us split the outgoing microwave into two equal parts 
and use phase-sensitive amplifiers with ip = Q and ip = 
7r/2 in the two channels. To keep the same noise S per 
channel, we increase the gain by the factor \/2, that also 
compensates the signal loss at the splitter. Then the 
channel ip — Q produces the "spooky" back-action (pS)) . 
while the channel ip = 1^/2 informs us of the "realistic" 
back-action (P^)) . and the relation ([50)1 between K and 



AI is the same as between K^ax and Al^ax in (1211) . 

In the ideal case (r/ = 1) the qubit evolution description 
can be translated into the language of the POVM-type 
measurement. In the same way as in Sec. |lll Eqs. ((28)) 
and (|29l) can be converted into the measurement operator 



M{Im,Qm) = CXp{-iKQmtcrj2) 



P\O){Im)\0){0\ + X^P\l)iIm)\l){l\ 



(35) 



where now S is the total (experimental) noise per quadra- 
ture, AI is the experimental response, and F = 8x'^v/k-\- 



where the probabilities Py-^ are given by Eq. ([3]). Simi- 
larly, Eqs. (US]) and (fTOj) for the case of a phase-sensitive 
amplifier can be converted into the same measurement 
operator (j35p . in which Qm is replaced with /„. 
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V. CONCLUSION 

We have presented a simple physical picture of the 
qubit evolution due to its measurement in the circuit 
QED setup. The "spooky" back-action is universal, it 
is caused by gradual extraction of information about the 
qubit state. The "realistic" back-action is due to a spe- 
cific mechanism: fluctuation of the photon number in 
the resonator. For a phase-sensitive amplifler the qubit 
evolution is described by Eqs. ^TE\\ and ([25)) : it is deter- 
mined by the output signal I{t), which corresponds to 
the amplified quadrature. For a phase-preserving ampli- 
fier the evolution is described by Eqs. (1251) and p2l) : it is 
determined by two output signals: I{t) and Q{t), where 
I{t) now corresponds to the quadrature, which provides 
information about the qubit state (S-quadrature in the 
resonator) and Q{t) corresponds to the orthogonal A- 
quadrature, which gives us a record of the photon number 
fluctuations in the resonator. 

While the cQED setup significantly differs from both 
the broadband quantum measurement setup psj and the 
standard optical setup [ill i see that the description of 
the qubit evolution is exactly the same as in both these 
cases if a phase-sensitive amplifier is used. The descrip- 
tion is only slightly different when a phase-preserving am- 
plifier is used: we should assign the "spooky" and "re- 
alistic" back-actions to the separate output signals I{t) 
and Q{t) instead of only one signal. It is also useful to 
think about the phase-preserving case via the model in 
which we split the outgoing microwave (the quantum sig- 
nal) into two equal parts and then use 90-degree-shifted 
phase-sensitive amplifiers for these two channels. 

We intentionally considered only the simplest case, be- 
cause most of the further steps and generalizations are 
quite straightforward 23]. In particular, it it very simple 
to include Rabi oscillations and energy relaxation of the 
qubit state. For that we have to take time-derivative of 
the evolution equations and add the terms due to Rabi 
oscillations and energy relaxation. If the Stratonovich 
definition of the derivative is used, we get the equations of 
the Bayesian formalism [23']; if the Ito derivative is used, 
we get the equations of the quantum trajectory formalism 
[1, [3| • Generalization to measurement of several entan- 
gled qubits is also straightforward ^2^. We have consid- 
ered only the resonant case, a;„i — uJr', however, gener- 
alization to the case ujm 7^ is quite simple (see [^): 
we just need a different definition of the informational B- 
quadrature and photon-fluctuation ^-quadrature. In our 
formalism we implicitly assumed sufficiently wide band- 
width of the amplifier (much larger than the ensemble 
dephasing F and Rabi frequency ^Ir). If this is not the 
case, the formalism should change significantly. How- 
ever, we believe that in most of the practical cases we 
can take this effect into account by adding a classical 
narrowband filter to the classical signal at the amplifier 
output; this will correspond to passing the signals I(t) 
and Q(t) through the low-pass filters. A much more seri- 
ous change of the theory is required when the resonator 



bandwidth n is comparable to F or flu; this still has to 
be done. 

Understanding the difference between the "spooky" 
and "realistic" back-actions is important for designing 
the quantum feedback control of the Rabi oscillations 
[r3 |. The simplest case is when a phase-sensitive am- 
plifier amplifies the informational i3-quadrature. Then 
there is no "realistic" back-action, and the feedback loop 
should only modulate the amplitude of the Rabi drive 
(i.e. the Rabi frequency ^R}i this case was well studied 
for the broadband setup [12|. The situation is differ- 
ent for a phase-preserving amplifier. Then we need two 
feedback channels: the first (usual) channel should mod- 
ulate the Rabi frequency iln to compensate the "spooky" 
back-action, while the second channel should compensate 
the "realistic" back-action by modulating the qubit fre- 
quency ujqh or the the frequency of the Rabi drive ujr. 
The controller for the second feedback channel is quite 
simple: it should compensate the contribution iKQ{t) to 
the qubit phase derivative (t>{t) due to the K-terva in Eq. 
([52)) . Therefore the controller is 

A(c^,b-c^j^) = -i^[Q(t)-(Q)], (36) 

i.e. we should directly apply the signal Q{t) to modulate 
LOqb or bjji. The second feedback channel essentially elimi- 
nates the K-leum in Eq. ()32p and decreases the ensemble 
dephasing F by K^S/i ^ (A/)^/4S'. Correspondingly, 
in the absence of the first (main) feedback channel the 
peak-to-pedestal ratio of the Rabi peak increases from 
2 to 4 in the quantum-limited case. The first feedback 
channel should be the same as for the broadband setup; it 
depends on the signal I{t) and can be realized using vari- 
ous ideas for the controller ( "direct" , Bayesian, "simple" , 
etc. [T^V Notice that without the second channel the 
feedback performance is determined by the quantum effi- 
ciency f), while with the second channel it is determined 
by 7) = i)/ {l — rj + f}) (this is one more combination of the 
terms in Eq. ^ , which can be used for the definition [2^ 
of quantum efficiency). The case of a phase-sensitive am- 
plifler, which amplifies a non-optimal quadrature {(p ^ 0, 
(f 7^ 7r/2) is similar to the case of a phase-preserving am- 
plifier, but both feedback channels should start with the 
same signal I{t). In both the phase- sensitive and phase- 
preserving setups a perfect feedback control is possible in 
the quantum-limited case rj — 1. 

Discussion of the "spooky" and "realistic" back-actions 
in the cQED setup necessarily raises the question of 
causality. When the microwave leaves the resonator, 
it does not yet "know" in which way it will be mea- 
sured (phase-preserving or phase-sensitive, which angle 
if, etc.). Moreover, when a circulator is used for the 
outgoing microwave, the field in the resonator and the 
qubit can never "know" in a realistic way which method 
of measurement is used. Nevertheless, the qubit evolu- 
tion strongly depends on the measurement method. As 
we discussed, the "spooky" evolution moves the qubit 
state along the meridians of the Bloch sphere, the "real- 
istic" back-action moves the state along the parallels, and 
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the measurement method determines whether the qubit 
experiences the "spooky" or "reahstic" back-action (or 
their combination). In this sense the "reahstic" back- 
action is not fully realistic: it has the physical mecha- 
nism, but whether this mechanism works or not is de- 
termined in a spooky way. The causality requires that 
we cannot pass a "useful" information to the qubit by 
choosing the measurement method. This means that the 
ensemble-averaged evolution of the qubit cannot depend 



on the measurement method (this is the general require- 
ment of causality in quantum mechanics). It is surely 
satisfied in our cQED setup. 
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